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A simple five-dimensional brane world model is proposed, motivated by M-theory compactified 
on a six-dimensional manifold of small radius and an S 1 jZ-x of large radius. We include a leading- 
order higher curvature correction to the tree-level bulk action since in brane world scenarios the 
curvature scale in the bulk may be comparable to the five-dimensional Planck scale and, thus, 
higher curvature corrections may become important. As a tractable model of the bulk theory we 
consider pure gravity including a (Ricci-scalar) 4 -correction to the Einstein-Hilbert action. In this 
model theory, after a conformal transformation to the Einstein frame, we numerically obtain static 
solutions, each of which consists of a positive tension brane and a negative tension brane. For these 
solutions, we obtain two relations between the warp factor and the brane tensions. The existence of 
these relations implies that, contrary to the original Randall-Sundrum model, the so called radion 
is no longer a zero mode. We conclude that the tension of our brane should be negative and that 
fine-tuning of the tension of both branes is necessary for a large warp factor to explain the large 
hierarchy between the Planck scale and the electroweak scale. 



I. INTRODUCTION 



The idea that our world may be a brane embedded in a higher dimensional spacetime has been attracting a great deal 
of physical interest. This idea is often called the brane world scenario and, as suggested by Randall and Sundrum 
may be able to explain the large hierarchy between the Planck scale and the electroweak scale in a natural way. Many 
aspects of the brane world scenario have been investigated: for example, the effective four-dimensional Einstein's 
equation on a positive tension brane ||, weak gravity j|]5), black holes ||, inflating branes 0, cosmologies [p|-^2|, 
and so on. 

In the original Randall-Sundrum brane world scenario, these authors considered five-dimensional pure gravity 
described by the Einstein-Hilbert action with a negative cosmological constant and two 3-branes with tension. Because 
of the negative cosmological constant the five-dimensional bulk geometry is highly curved and the curvature scale 
is possibly of the order of the five-dimensional Planck scale, while the induced geometries on the branes are flat, 
provided that brane tensions are fine-tuned. Therefore, it is expected that quantum effects in the bulk may be 
important in the brane world scenario, since quantum effects in curved spacetime usually become important when 
the spacetime geometry is highly curved or when the causal structure is non-trivial |13|. In this connection, several 
authors investigated quantum effects in the brane world scenario [T^-[T^]. 

In particular, in ref. |l6[ ] exact semiclassical solutions representing a static brane world with two branes were 
obtained by analyzing the semiclassical Einstein's equation in five-dimensions with a negative cosmological constant 
and conformally invariant bulk matter fields. There, the following two types of solutions were found. Type-(a): 
solution with a positive tension brane and a negative tension brane. Type-(b): solution with two positive tension 
branes. For each type of semiclassical solution, two relations between the warp factor and brane tensions were found: 
one giving the warp factor as a function of the brane tensions and another giving a relation between the brane tensions. 

Although it is interesting that we could obtain analytic solutions in the model of ref. , it seems that this model 
is not realistic enough. As far as the author knows, there is no realization of conformally invariant bulk matter fields 
starting from M-theory or superstring theory. Nonetheless, it is expected that the solutions in ref. [jl6f may actually 
capture some important features of quantum effects in the brane world scenario. Hence, it is worth while to extend 
the analysis of ref. jl6| to more realistic brane world models which also take bulk quantum effects into account. For 
this purpose, one would like to consider higher curvature corrections to the tree-level bulk action. As discussed in the 
next section, i? 4 corrections are realistic from the point view of M-theory. 

In this paper a simple five-dimensional brane world model is proposed, motivated by M-theory compactified on a 
six-dimensional manifold of small radius and an S 1 jZ-i of large radius. We include the leading-order higher curvature 
correction to the tree-level bulk action. As a tractable model of the bulk theory we consider pure gravity including 
a (Ricci-scalar) 4 -correction to the Einstein-Hilbert action. In this model theory, after a conformal transformation 
to the Einstein frame, we numerically obtain static solutions, each of which consists of a positive tension brane and 
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a negative tension brane. For these solutions, we obtain two relations between the warp factor and brane tensions. 
Those solutions and relations are a close analogue of the type- (a) solutions and relations obtained in the model of 
ref. On the other hand, in the present model it will be shown that there is no analogue of the type-(b) solutions. 
This fact might be considered to be consistent with the suggestion of refs. [|l4|Jl8| that the type-(b) solutions are 
unstable. 

This paper is organized as follows. In Sec. || we describe a simple brane world model which take bulk quantum 
effects into account. In Sec. Ill we numerically obtain static solutions in the model, and two relations between the 
warp factor and brane tensions are derived. Sec. |^ is devoted to a summary of this paper. 



II. MODEL DESCRIPTION 



In this section we propose a simple brane world model, motivated by M-theory compactified on a six-dimensional 
compact manifold (eg. Calabi-Yau manifold) of small radius and an S 1 jZ<x of large radius In this situation, 

effectively we may consider a five-dimensional theory compactified on the S 1 jZi. In the five-dimensional bulk action, 
we shall consider a correction by a R 4 term to the Einstein-Hilbert term [] since several calculations of higher-order 
corrections to the effective action suggest that in eleven-dimensions R 2 terms do not appear but R 4 terms may 
appear pl],^2| 0. It is expected that the R 4 corrections play important roles in brane world scenarios since the 
curvature scale in the bulk may be comparable to the five-dimensional Planck scale and, thus, higher curvature 
corrections cannot be neglected. 

On the other hand, as for the action on branes, higher curvature corrections are expected to be less important than 
those in the bulk and can be neglected since curvature scale induced on branes (at least on our brane) should be small 
compared to the Planck scale in low energy. Nonetheless, motivated by the four- and ten-dimensional effective theory 
induced on the fixed points of S 1 /Z2 |U9lE4|, we may include R 2 corrections to brane actions. In the following, we 
will explicitly see that the R 2 corrections do not play any roles in our analysis of static solutions. 

Since general higher curvature terms are difficult to analyze, for simplicity, we shall consider Ricci scalars only. 
Furthermore, we assume that the compactification from eleven dimensions to five dimensions is properly stabilized 
and, for simplicity again, we do not consider the corresponding moduli as dynamical fields in five dimensions. Namely, 
in our analysis, we shall consider the following action. 



2k- 



: R 5 + clk 2 R 4 



A 



d 4 y^q(bR 2 4 — A) + / d^y/^qibR 2 - A) 



(1) 



where k and A are the five-dimensional gravitational constant and cosmological constant, a, b and b are dimensionless 
constants, and A and A are brane tensions. The fixed-point hypersurface, or the world volume of a 3-brane, S (or £) 
is represented by x 
defined by 



Z M (y) (or x — Z M (y), respectively) and the induced metric q^ v (or g M „, respectively) is 



%v{v) = e M {y)e v (y)gMN\ x=Z (y), 
8Z M (y) 



GO 



(2) 



(or q^(y) = eff (y)e£ (y)g M N\ x =z(y)i e fl (y) = dZ M (y)/dyi*, respectively). The Ricci scalars i? 5 , R A and R± are of 
Qm n , q^v and q^, respectively. 

Following ref. pa], we perform the conformal transformation 



9m N 



e> ' ' yMN , 



Kip = —= ln(l + 8aK 4 R 3 5 ) 
v3 



(3) 



to obtain the following expression. 



1 We, off course, include the Einstein-Hilbert term since it appears in the tree- level effective action in eleven dimensions [20|. 
The importance of R 4 terms in M-theory was originally pointed out in ref. and many authors showed evidences of it [22|. 



Possible cosmological consequences were discussed in ref. 
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/ = / d 5 x 



M 



d yy/-q 



L 



d y\ -q 



^-\g MN d M ^d N ^-u^) 



b (r 4 + V3KD 2 ip - ^q^d^dvip 



b ( i?4 + V3kD tp — \^?q^ V dfj,l/)d v lp 



(4) 



where the conformally transformed induced metrics q^ v and q„ v are defined by 



%M = eU^VWq^y) = e™(y)e»(y)g M N\ x =z(y) 
LM = eW^&Wq^y) = ef {y)e? {y)g MN \ x = 



(5) 



D and D are four-dimensional covariant derivatives compatible with <7 M „ and q^, respectively, and the Ricci scalars 

i?5, i?4 and i?4 are of §mn, q^v and g^,,, respectively. The potential U(ip) and the functions f{ip) and /(VO are given 
by 



A + (3/16) K ~ 10 / 3 a^ 3 (e ( ^ 2 ^ - l) 4 / 3 

f(1>) = e -(2/V3)KV- A) 

/(V) = e ~( 2 /^)^A. (6) 
To obtain the expression (^) we have assumed that 8k 4 aR^ > — 1. For negative A, the potential can be rewritten as 

Uty) = \A\e-^ 6 ^ [-1 + a(e^ 2 ^ - l) 4 / 3 ] , (7) 

where 

a=(3/16)^ 10 / 3 a- 1 / 3 |A|- 1 . (8) 

Since the ^-dependent part of the action (Q) is of the Einstein-Hilbert form, the conformal frame in which the metric 
is gMN is called Einstein frame. On the other hand, we shall call another conformal frame in which the metric is qmn 
the original frame. 

In this paper, we assume that A < and consider a static configuration with the ansatz 



g MN dx M dx N = e- 2A ^f 1 ^dx' 1 dx u + dw 2 , 
ip = ip(w). 



(9) 



This ansatz represents a general configuration with the four-dimensional Poincare invariance. Off course, the set of 
all configurations with the four-dimensional Poincare invariance in the Einstein frame is equivalent to that in the 
original frame. With this ansatz the curvature-squared term in the brane action does not contribute to the equation 
of motion at all. Einstein's equation and the field equation of the field ip are given by 



d 2 A 
dw 2 



dw 



!r> = 11 - 

dw 



„4A 



dw 



dw 
dw 



U'(tp) = 0. 



(10) 



Note that the third equation is dependent of the first two equations unless dip/dw = (the Bianchi identity), while 
the first equation can also result from the second and the last equations unless dA/dw = 0. When we compactify the 
w-direction by S 1 /Z? so that w ~ w + 2L and that w ~ — W, there appears the following matching condition for ip. 
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u>— >Z— dill iw—Ij 

where we have assumed that £ and £ are world- volumes of branes at the two fixed points w = and w = L, 
respectively. We suppose that the brane at w = is our world and shall call it our brant. We shall call another brane 
at w — L the hidden brane. As for the function A(w) in the metric, we have the junction condition 

dA 

6 lim — = K 2 f(ip)\ m 
dA 

6 lim =- K if W \ (12) 

This is a special case of Israel's junction condition j2(|. In the Einstein frame the so called warp factor can be defined 
by 4>e = e A (°^/e A ( L \ Correspondingly, the warp factor in the original frame is 



= exp \ [A(0) - A(L)} + ^= MO) - tl>(L)] }. (13) 



It is evident that ip = ipQ is not a solution because of the matching condition ( |11| ) , where ipo is an extremum of the 
potential U(tp). In particular, we can show that 

lim U(ip) = lim UU) = 0, (14) 

w— >L— 



and thus r/> cannot stay at tpo unless A is zero. Actually, provided that equations of motion ( |10[) and the junction con- 
dition (|l^) are satisfied, the matching condition (|ll|) is equivalent to the vanishing-potential condition ( |l4] ) combined 
with 

dA dib 
lim r- < 0, 

in— >+o dw 

lim r- < 0. (15) 

w—>L-o dw dw 



In order to see the necessity of the condition (|15[), notice that /(V0/'(V0 < and / (V>)/'(V0 < 0. 



III. NUMERICAL SOLUTION 



For the purpose of numerical integration, it is convenient to rewrite all equations in terms of dimensionless variables. 
Hence, we introduce the dimensionless independent variable x defined by x — w/L and consider the region < x < 1, 
where L is the distance between two branes. As for the dependent variables, we introduce the following three: 



2/1 0) 

J/2 (a:) 
2/3 (z) 



A, 

L dA 

dw ' 

K 

2V3 



Differential equations for these dimensionless independent variables are given by 



(16) 



2/1 = 2/2, 

y 2 =4[yi + (L/l) 2 V(y 3 )], 
(h) 2 = Vl + (VOMlte), 



(17) 



where dots denote differentiation with respect to x, the length scale I is defined by I = k 1 y6/|A|, and 



v% 3 



-l + a(e 33 ' 3 - I) 4 / 3 



(18) 
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As already mentioned in the previous section we assume that A < 0. The set of these three differential equations is 
equivalent to the equation of motion (|l^) as long as 1/2 is not zero. As we shall argue later, there is no static solution 
for a < 1. Hence, for the present we shall concentrate on the case a > 1 only. The potential is shown in Figure |l| for 
a = 2.0, 1.5, and 1.2. The vanishing-potential condition ( |l4|) is written as 

V(y 3 (0)) = V(y 3 (l)) = 0, (19) 

and should be complemented by 

y 2 (0)y 3 (0) < 0, 

y 2 (i)y 3 (i) < 0. (20) 

It is easy to impose the boundary condition ( |l9| ) since the roots of V(jj 3 ) are analytically obtained as 2/3 = y± for 
a > 1, where 

y ± = iln(l±a- 3 / 4 ). (21) 

The complementary condition (^0|) should be checked after a solution of the differential equation ([11]) with the 
boundary condition (n9) is obtained. Thence, the junction condition (n2[) determines the brane tensions as 



X/(6K- 2 r 1 ) = (l/L) y2 (0)e 4 y^°\ 

X/iGn- 2 !- 1 ) = -(l/L)y 2 (l)e 4 y^\ (22) 
Finally, the warp factor <f> given by ( |l3| ) is written as 

<t> = exp M0) + 2, 3 (0) - y x (l) - 2/3(1)] . (23) 
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Note that, without loss of generality, we can impose the additional condition 



2/1(0) = o, 



(24) 



since none of the above equations is changed by the shift y\ (x) — > y\ (x) 
with the vanishing-potential condition (n9J) , which can be rewritten as 



2/3(0) 

2/3(1) 



V±, 

y±, 



yi (0). This additional condition combined 



(25) 



give enough number of boundary conditions for the set of three differential equations ([17]). Here, plus or minus signs 
in two of (|2^) are independent. According to four possible choices of the signs in (25), there are four possible types 
of solutions. 



(++)-type 
(+-)-type 
(-+)-type 

( — )-type 



2/3(0) 


= y+ 


2/3(1) 


= 2/+, 


2/ 3 (0) 


= y+ 


2/3(1) 


= y-, 


2/3(0) 


= y- 


2/3(1) 


= y+, 


2/3(0) 


= y- 


2/3(1) 


= y-- 



(26) 



We shall solve the differential equation (|T^) with the boundary condition (^4|) and ( |25| ) by the so called relaxation 
method j2^] . This method works very well if a good initial guess is given. In the following we shall solve the differential 
equation many times, each time with different values of the parameters a and L/l. Hence, the previous solution can 
be used as a good initial guess in the next calculation with slightly different parameters. 

Figures pi ||, [|, ||, || and ^| show several ( — h)-type solutions obtained by the relaxation method. Physical parameters 
in these solutions are a = 2.0, 1.5, 1.2 and L/l = 10.0, 20.0, 30.0. 




FIG. 2. The numerical solution yi (x) for L/l = 10.0 and a = 2.0, 1.5, 1.2. 
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L/l=20.0 




FIG. 3. The numerical solution yi (x) for L/l = 20.0 and a = 2.0, 1.5, 1.2. 



L/l=30.0 
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x 

FIG. 4. The numerical solution y t (x) for L/l = 30.0 and a = 2.0, 1.5, 1.2. 
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FIG. 5. The numerical solution y 3 (x) for L/l = 10.0 and a = 2.0, 1.5, 1.2. 
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FIG. 6. The numerical solution y 3 (x) for L/l = 20.0 and a = 2.0, 1.5, 1.2. 
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FIG. 7. The numerical solution y 3 {x) for L/l = 30.0 and a = 2.0, 1.5, 1.2. 



From these figures, we can see that, except for vicinities of the boundaries, j/3 stays near the minimum y n 
potential V(ys) and y\ is almost linear in x: [] 



L 



yi(x) ~ - — \J\V(y m in)\ x, 



of the 



(27) 



In other words, except for vicinities of the boundaries, the scalar field ip stays near the minimum of the potential U(ip) 
and the five-dimensional geometry is almost the AdS whose curvature is determined by the minimal value of U^). 
However, because of the boundary condition ((2^), near-boundary behaviors of solutions are non-trivial. Actually, 
figures H and ^ show that the five dimensional geometry near the boundaries deviates rather strongly from AdS. The 
deviation is larger for a smaller value of a, or a larger value of the coefficient a of the i? 4 -term, as easily expected. 
Note that y^ — y\ would be independent of x if the geometry was AdS. Each of figures and 11 shows how the 
scalar field approaches one of roots of the potential. 



3 As we shall argue in the paragraph after the next, these approximate expressions could be inferred without any numerical 
calculations for large values of L/l. 
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L/l=1 0.0 
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FIG. 8. The numerical solution 2/2(1) = yi(x) in the vicinity of the boundary x — for L/l = 10.0 and a = 2.0, 1.5, 1.2. 



L/l=10.0 



CM 




FIG. 9. The numerical solution 2/2(2:) = yi(x) in the vicinity of the boundary x = 1 for L/l = 10.0 and a = 2.0, 1.5, 1.2. 
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FIG. 10. The numerical solution 2/3(2;) in the vicinity of the boundary x = for L/l — 10.0 and a = 2.0, 1.5, 1.2. 
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FIG. 11. The numerical solution ys{x) in the vicinity of the boundary x = 1 for L/l — 10.0 and a = 2.0, 1.5, 1.2. 

We now argue that there are no (static) solutions of (++)- and ( )-types. First, let us rewrite the differential 

equation as 

2/i = V2, 
2/2 = iyl, 



where 



2/3 = 2/4, 

2/4 = 4 1 / 2 y 4 + (i/0 2 ^'(2/3)/2, (28) 



2/4^) = — (29) 



and the corresponding boundary condition is 



It is easy to show from (B8f) and (pG) that 



2\/3 dw' 



2/i (0) =0, 
2/ 2 (0) + j/ 4 (0) =0, 
2/3(0) = y ± , 

2/2(1) + 2/4(1) =0. (30) 



1 

dxe- 4va V'(y 3 ) = 0. (31) 





Next, it is also easy to show by using the last equation of ( ps| ) that, if yi{x\) < and 2/3(^1) < j/mm for < 3 xi < 1, 
then 2/4(2;) < for x\ < v x < 1, where j/mm is the global minimum of ^(2/3) between y_ and 1/+. Thus, if 2/3 starts from 
2/+ at 2; = and reaches the region 2/3 < y m in then 2/3 cannot return to 2/+ at x = 1 . Thirdly, combining this fact with 
(31), we can show that there is no solution of the (H — h)-type which is bounded in the region y~ < ys < y+. Actually, 
(31) implies that if 2/3 starts from y + at x — and is bounded in the region y~ < y% < y+ then 2/3 should reach the 
region 2/3 < y m in > from which 2/3 cannot return to y+ at x = 1 . Similarly, we can show that there is no solution of the 
( )-type which is bounded in the region 2/- < 2/3 ^ 2/+- Heuristically, the above mathematical statement can be 



easily inferred if one interprets the last two equations of (28) as the equation of motion for a particle position 2/3(2;) 
at time x. In this interpretation, the particle receives the force due to the reversed potential —{L/l) 2 V{y^)/2 and 
the friction (or anti-friction) force 4?/22/3- Moreover, with this heuristic interpretation in mind, it seems likely that 
any solution should be bounded in the region y~ < yz < y+ since the reversed potential — {L/l) 2 V{yz)/2 is negative 
outside the region j/- < 2/3 < y+ and V(ys(0)) — ^(2/3(1)) = 0. Finally, combining this with the above mathematical 

statement, we can say at least for moderate a that there are no solutions of (++)- and ( )-types. 

Moreover, (H — )- and ( — (-)-types are equivalent to each other up to the coordinate transformation x — > 1 — x. 
Hence, (— +)-type solutions are all we have been seeking. 



Figures [L2| and |13| show relations between the warp factor <f> given by (23) and the brane tension A or A given by 
( p2|) for ( — h)-type solutions. We can see that the tension of each brane converges quickly to an a-dependent value 
when the warp factor becomes large. To be precise, 

-0.54 (a = 2.0) 

A/(6K" 2 r 1 ) -> { -0.36 (a = 1.5) , (32) 
-0.18 (a = 1.2) 

1.28 (a = 2.0) 

A/(6/t" 2 r 1 ) -> { 1.34 (a = 1.5) (33) 
1.40 (a = 1.2) 

as (f> — > 00. The convergence can be easily understood as follows by the heuristic interpretation. We consider a 
particle in one dimension whose position at time x is 2/3(2;) and which receives a force due to the reversed potential 
— (L / l) 2 V (yz) / 2 and the friction (or anti-friction) force. The particle moves from 2/3 = 2/- 1° 2/3 = V+ m the fixed 
duration < x < 1. For a large L/l, in order to satisfy this boundary condition, the particle should stay near 2/3 = y m in 
for a relatively long time since the reversed potential is steep. In this case, the initial (or final) velocity 2/4(0) (or 
2/4(1), respectively) should be fine-tuned to a value close to the 'escape velocity', which is roughly proportional to 
L/l, against the reversed potential. Because of the relation 2/2(0) + 2/4(0) = (or 2/2(1) + 2/4(1) = 0, respectively) 
and the junction condition (|22|), the fine-tuning of 2/4(0) (or 2/4(1): respectively) is equivalent to the fine-tuning of 
A (or A, respectively). The required value of A (or A, respectively) is almost independent of L/l since the required 
value of 2/4(0) (or 2/4(1)1 respectively) is roughly proportional to L/l as stated above. On the other hand, when 2/3 
stays near i/ mm with a very small velocity, 2/2 should satisfy y 2 ~ —(L/l) 2 V(y m i n ) since the right hand side of the 
last equation of (|l7|) should vanish approximately. Hence, 2/1 grows approximately linearly in x with the growth rate 
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proportional to L/l when 1/3 stays near y m in- Thus, the exponent of the warp factor should be roughly proportional 
to L/l. Actually, except for vicinities of the boundaries, the numerical solutions are well approximated by ( p7|) and 
thus the warp factor is well approximated by 



exp 



jV\ v (ymin)\ +y- -y+ 



(34) 



Finally, combining the approximate linearity of the exponent of the warp factor with the fine-tuning of the brane 
tension, we can conclude that the brane tension should converge to a constant as the warp factor becomes large. The 
limiting value of the brane tension should depend on the parameter a since the 'escape velocity' depends on a. This 
conclusion is of course consistent with the numerical result. 

Figure [l4| shows a relation between tension of our brane at x = and that of the hidden brane at x = 1. This 
relation can be considered as a necessary condition for the system with two branes to be static, or the condition for 
the four-dimensional cosmological constant to vanish. 
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FIG. 12. The relation between the warp factor <f> given by (|23|) and the tension A of our brane at x = 0, which is given 
by (p^), for the ( — +)-type solutions. The horizontal axis represents ln<^ and the vertical axis represents A/(6k -2 / -1 ). The 
physical parameter is a = 2.0, 1.5, 1.2. As cj> becomes large, A/(6k~ 2 Z -1 ) converges quickly to the a-dependent value given by 
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FIG. 13. The relation between the warp factor (f> given by (J23j) and the tension A of the hidden brane at x = L, which is 
given by (|22|, for the ( — +)-type solutions. The horizontal axis represents ln0 and the vertical axis represents A/(6fi; _2 £ _1 ). 
The physical parameter is a = 2.0, 1.5, 1.2. As <j> becomes large, \/(6k~ 2 1~ 1 ) converges quickly to the a-dependent value given 
by ph. 



Constraint between Tensions 




Dimensionless Tension of Hidden Brane 



FIG. 14. The relation between tension A of our brane at x = and the tension A of the hidden brane at x = 1. The 
horizontal axis represents A/(6k -2 / -1 ) and the vertical axis represents \/(6k~ 2 1~ 1 ). This relation can be considered as a 
necessary condition for the system with two branes to be static, or the condition for the four-dimensional cosmological constant 
to vanish. The physical parameter is a = 2.0, 1.5, 1.2. 
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Now let us give an argument that there is no static solution for a < 1. The following argument is almost the same 

as that for the non-existence of (++)- and ( )-type solutions. First, let us consider the case < a < 1. In this case 

the potential 1^(2/3) has only one root y = y + , and thus only (+- |-)-type solutions are allowed if any. However, since 
V(ll3) > for 2/3 < y + , the condition (|3l| ) excludes (++)-type solutions which are bounded in the region 2/3 < y+- 
On the other hand, since the reversed potential — (L / l) 2 V (y^) / 2 is negative for 2/3 > y + and V{yz(Q)) — V(yz{l)) = 0, 
the above heuristic interpretation in terms of a particle motion strongly suggests that, if there exists a solution of the 
differential equation, the solution should be bounded in the region 2/3 < ?/+■ Hence, we can conclude that there is 
no solution for < a < 1. Next, for a < 0, there is no root of the potential 1^(2/3). Hence, in this case there is no 
way to satisfy the boundary condition (|25|). Therefore, there is no solution for a < 1, and it is actually enough to 
concentrate on the case a > 1 as we did. 



IV. SUMMARY AND DISCUSSIONS 



We have proposed a simple five-dimensional brane world model, motivated by M-theory compactified on a six- 
dimensional manifold of small radius and an S 1 /Z 2 of large radius. We have included the leading-order higher 
curvature correction to the tree-level bulk action since in brane world scenarios the curvature scale in the bulk may 
be comparable to the five-dimensional Planck scale and, thus, higher curvature corrections may become important. 
As a manageable model of the bulk theory we have considered pure gravity including a (Ricci-scalar) 4 -correction to 
the Einstein-Hilbert action. 

In this model theory, after a conformal transformation to the Einstein frame, we have numerically obtained static 
solutions, each of which consists of a positive tension brane and a negative tension brane. The solutions are parame- 
terized by a dimensionless parameter a in the bulk theory and L/l, where L is distance between two branes and / is 
a length scale determined by the (negative) bulk cosmological constant. Several solutions are shown in Figures ^, ||, 
UJandg. 

The warp factor and tension of both branes have been calculated for various values of a and L/l and, by eliminating 
L/l, we have obtained two a-dependent relations between the warp factor and brane tensions. The existence of these 
relations implies that, contrary to the original Randall-Sundrum model, the so called radion is no longer a zero mode. 
In this sense, the present model is similar to those in Refs. |2^j2^||. The two relations completely determine the 
brane tensions as functions of the warp factor and are shown in Figures ^ and From these figures we conclude 
that the tension of our brane should be negative and that fine-tuning of the tension of both branes is necessary for 
a large warp factor to explain the large hierarchy between the Planck scale and the electroweak scale. To be precise, 
the brane tensions should be fine-tuned with high accurac y t o values shown in equations ([32]) and (|33|). 

Further, eliminating the warp factor from Figures |l| and]l|, we have obtained a relation between the brane tensions. 
It is shown in Figure 14 and can be considered as a necessary condition for the system with two branes to be static, 
or the condition for the four-dimensional cosmological constant to vanish. Namely, unless this relation is satisfied, the 
system cannot be static but becomes dynamical, regardless of initial conditions (i.e. initial position of branes, initial 
velocity of branes, and so on). 

A stability analysis of solutions obtained in the present paper is an important topic for future work. Here, we 
only offer a comment concerning this subject: we cannot derive a correct effective action by simply substituting the 
solutions into the action. Actually, if we substitute any static solutions into the action then the action vanishes 
because of the Hamiltonian constraint. A simple illustration of this fact is given in Appendix [A|. 

Several extensions of the present model may also be of interest for future work: (i) inflating brane solution; (ii) 
cosmological solution; (iii) inclusion of Ricci tensor and Weyl tensor contributions to the i? 4 -tcrm in the action; (iv) 
inclusion of a 3-form field and modulus corresponding to the six-dimensional compactification. (i) It is probably not 
difficult to extend the static solutions in the present paper to inflating brane solutions. The relation between brane 
tensions, corresponding to Figure |lj, is expected to become dependent on the four-dimensional cosmological constant 
induced on branes as well as the model parameter a. (ii) Extension to cosmological solutions should be possible. This 
should not be as difficult as extension of the semiclassical solutions in ref. jl6[ to cosmological solutions. The latter 
seems rather difficult because of the so called moving mirror effect EJ, which is non-local. On the other hand, the 
present model has a locally defined Lagrangian density in five dimensions. Hence, extension to the cosmological context 
is easier in the present model than in the model of ref. |lqj . Moreover, the present model seems more realistic and, thus, 
worth while investigating in more detail, (iii) Although we have investigated effects of the (Ricci-scalar) 4 -correction 
only, it would be desirable to investigate effects of other forth-order curvature terms. Note, however, that effects of 
forth-order Weyl terms are probably less important insofar as we consider the metric (^|) or small perturbations around 
it, since the Weyl tensor vanishes for the metric ([)]). (iv) In the present model we have considered pure gravity in the 
bulk. However, more realistic model should include a 3-form field existing in the bosonic sector of eleven dimensional 
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supergravity as well as moduli fields due to the compactification from eleven dimensions to five dimensions. In this 
case, we may consider effects due to various eighth-order derivative terms other than i? 4 terms. 
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manuscript. He would be grateful to Professor J. Yokoyama and Dr. T. Shiromizu for helpful discussions in the early 
stage of this work. In particular, the author would greatly appreciate Professor J. Yokoyama for informing as to 
his interesting work on topological i? 4 inflation |2^]. This work is supported by CITA National Fellowship and the 
NSERC operating research grant. 



APPENDIX A: ESTIMATE OF THE ACTION 



By the ansatz (||), the action (||) is reduced to 
/ / d 4 xC, 




cPA 5 (dA\' 
dw 2 \dw J 



where the integration with respect to w in this expression is over the whole S , and 

U(tl>) = U(i>) + f(ip)6(w) + m)S(w - L). 
From this reduced action, the following equations of motion are derived. 



(Al) 



(A2) 



dw 



d 2 A 
dw 2 



0. 



dw V dw 



U'(ifj) = 0. 



(A3) 



These are equivalent to equations (|T(]), ( |TT| ) and (|l^), provided that the identification w ~ w + 2L ~ — L is imposed. 
We can estimate the value of the reduced action by using these equations. Actually, the first of (A3) reduces the 
effective Lagrangian density C to 



d 2 A 
dw 2 



dA 

dw 



dw J 



(PA 
dw 2 



dw-^- | e 

dw 



-4A 



dA 
dw 



= 0. 



(A4) 



Therefore, the effective Lagrangian density vanishes if a solution of equations of motion is substituted. This fact can 
be easily understood as follows in the static case without boundaries the Lagrangian of the system is just minus 
the Hamiltonian, which should vanish because of the Hamiltonian constraint. 



The above arguments can be applied to the situation in ref. 28 by simply replacing U(iJ)), f(tp) and f(ip) with 
appropriate functions. Since the above analysis indicates a vanishing effective potential, it is impossible to obtain 
a correct effective potential for the so called radion, which corresponds to L in the above arguments, by simply 
substituting solutions to the action. It seems that the non- vanishing effective potential obtained in ref. |28| merely 
measures an amount of inconsistency in their analysis. Nonetheless, their conjecture that the radion can be stabilized 
by inclusion of a bulk scalar held seems correct if the backreaction of the scalar field to the geometry is sufficiently 
small [S. 



The author would thank Professor J. Yokoyama for pointing out this. 
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